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Abstract 

We consider systems of two and three electrons in a two-dimensional 
parabolic quantum dot. A magnetic field is applied perpendicularly 
to the electron plane of motion. We show that the energy levels corre- 
sponding to states with high angular momentum, J, and a low number 
of vibrational quanta may be systematically computed as power series 
in 1/| J|. These states are relevant in the high-S limit. 
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1 Introduction 



Recently, few-electron systems in nearly two-dimensional semiconductor struc- 
tures has been a subject of intense theoretical and experimental research [1]. 
In the present paper, we continue a programme, initiated in [2, 3], aimed 
at providing reliable analytical estimates for the energy levels of model two- 
dimensional systems on the basis of a simple physical picture for the states. 

The idea is to use the inverse of angular momentum as an expansion 
parameter in Schrodinger equation, thus obtaining nonperturbative series 
for the energy and the wave function. Variants of this method have been 
widely used in Atomic and Nuclear Physics [4, 5]. 

The static exciton in a magnetic field was considered in [2]. Energy 
levels were computed as a function of the mass ratio and the magnetic field. 
Comparison with estimates obtained from two-point Pade approximants [6] 
yielded excellent results at any magnetic field strength. 

In paper [3], a model of three electrons in a quantum dot with 1/r 2 repul- 
sion, considered previously in [7], was studied. The 1/| J|-method provided 
analytical estimates for the energy and a set of approximate quantum num- 
bers to label the states. We will make use below of many of the results 
obtained in [3]. 

Systems of two and three electrons in a two-dimensional parabolic quan- 
tum dot are to be studied in the present paper. A magnetic field is applied 
perpendicularly to the dot's plane. In this problem, there is only one rele- 
vant parameter, (3 = (E cou i/TiVt) l l % , where E cou \ is the characteristic Coulomb 
energy, and Q is expressed in terms of the dot and the cyclotron frequencies 
as £1 = \JcOq + ujI/A. (3 — > will be called the oscillator limit, and (3 — > oo - 
the Wigner limit, suggesting that we are dealing with a few-electron version 
of the Wigner solid. In typical quantum dot experiments, (3 ~ 1. 

The physical picture emerging from the 1/| J|-expansion is the following. 
In the leading approximation, the electrons are located at the vertices of a 
regular polygon, which side minimises the total classical energy (centrifugal 
plus potential). The first corrections take account of small (harmonic) vi- 
brations around the equilibrium configuration. Higher corrections come from 
anharmonic oscillations. This picture is supported by numerical calculations 
for few-electron systems not only in the strong coupling regime [8], but in 
the region (3 ~ 1 as well [9]. 

A special comment deserves paper [9], in which the few-electron prob- 
lem in Eckardt coordinates is treated quasiclassically. Although the starting 
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points are far apart, our results are very similar up to the harmonic approx- 
imation. The main difference between paper [9] and ours is the expansion 
parameter, which is h in [9], and 1/|J| in our work. Higher order of the 
quasiclassical expansion are not reported in [9]. 

The plan of the paper is as follows. In section 2, we present the l/\J\- 
method. Results for two and three electrons are shown in sections 3 and 4. 
Concluding remarks are given in the last section. 

2 The 1/| J | -expansion 

We start with the dimensionless hamiltonian governing the two-dimensional 
motion of N electrons in a quantum dot of energy hu Q , and in a magnetic 
field normal to the plane of motion, 

The conventions follows, jj, is the electron effective mass, uj c = 

eB/fj,c is the cyclotron frequency, g is the effective gyromagnetic factor, S z 



is the z component of the total spin of the system, = JoJq + is the 



effective dot frequency. The length unit is yh/(fiil). The parameter (3 is 

given by f3 3 = \J -f^jp / Ml, where k denotes the dielectric constant. 

The actual parameter governing the problem is (3. When f3 « 1, the 
levels are small perturbations around oscillator energies. On the other hand, 
when f3 — > oo, the energy minimum is reached in a classical configuration 
which can be seen as a few-body model of a Wigner crystal. Both limiting 
situations conduce to exactly solvable problems. 

Through the introduction of Jacobi coordinates, the center-of-mass mo- 
tion is separated from the internal motion, 



H H crn Hi n t /<-.\ 



where 
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We have written ifj„ t in terms of rotational invariant coordinates, i.e., the 
moduli of the Jacobi vectors and the angles between them. The Jacobi 
vectors were defined as, 



Pk = vWpi jnt+i - £ £ 



n , fc=l,2,...,JV-l, (5) 



where the reduced masses are = k/(k + 1). The angle between and 
Pfc+i is denoted J labels the total internal angular momentum. 
The eigenvalues of H cm are simply, 

CYY~L i i C / \ 

~J~Q 2T?, cm + |t/cm| "I - 1 "I - ~^qJ cm- (6) 

We shall obtain approximate expressions for the eigenvalues of H int . We 
consider states with high \ J\ values, which are the relevant states at high 
magnetic fields. It is intuitively evident that at high \J\, E int ~ |J|, < p| >~ 
| J\. The solution of the Schrodinger equation may be organised as a power 
series in 1/|J|. A scaling of dimensions such that p 2 — > \ J\p 2 makes evident 
the dependence of each term on |J|. The scaled hamiltonian is written as 



h = — 



+ 



1 
J 

N-l 



H. 



VI it 



2n 



g^c q 
i 



d \ 



J 2 I \dp 2 k pkdpk, 
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V-:i | 



J k 

1 d 2 ) 2t ^ / 1 1 \ d 



N-2 I j 



\ 9^ 



(7) 
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where we have introduced the "renormalised" coupling constant (3 3 = (3 3 /\J\ 3 ^ 2 . 
When taking the formal limit |J| — > oo, f3 3 is kept fixed to take account of 
Coulomb repulsion nonperturbatively. 

In what follows, we consider only the two- and three-electron problems. 
The only term surviving in the r.h.s of (7) when |J| — > oo is the effective 
potential. Its absolute minimum gives a classical contribution to the energy. 
It is reached in a configuration where the particles sit at the corners of a 
regular polygon. In this configuration, the effective potential is a function of 
Pi only, 

2 sin 2 ir/N Npj gsin7r/7V 1 
e// Np\ + 8sin 2 vr/iV + Pl ^|sin%/2|' W 

where % is the angle between particles % and j, measured from the cm. Min- 
imisation of U e ff leads to an equilibrium value for pi, p w . The equilibrium 
values of the others coordinates, p k o, 9 k o, are determined from geometry. 

Higher contributions to the energy come from relaxing the equilibrium 
configuration, p k = p k0 + y k /\J\ 1/2 , 9 k = 9 k0 + z k /\ J\ 1/2 . The r.h.s. of (7) 
may, thus, be expanded as 

u u h 3 h A 

h = ho + jT\ W~ 2 W 2 " { ) 

where ho = U e ff(p w ), h 2 describes in general harmonic oscillations in ficti- 
tious magnetic fields, and h 3 , /i 4 , etc account for anharmonicities. 

For the spatial wave function and the scaled energy (related to h), we 
write series like (9) 



f = € » + |7| + ^ + --- ' 

These expressions are substituted into Schrodinger equation, leading to 
the chain of uncoupled equations, 



h = eo (12) 
h 2 ipo = e 2 ^ , etc. (13) 
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Higher corrections to the energy are obtained from ordinary perturbation 
theory, where h^, h 4 , etc. are interpreted as perturbations. 
Below, we present results for two and three electrons. 



3 Two electrons 

For two particles, there is only one Jacobi coordinate, p\ = r 2 — r\. It's 
equilibrium value (the modulus) satisfies the equation 

= "2 + \p\ Q - ^pio- (14) 

Writing p 1 = p w + y 1 /\ J] 1 / 2 , and substituting into the r.h.s. of (7), we 
obtain the operator coefficients I12, hs, etc. In particular 



d 2 I 

= u^ala,! + i), (15) 



where = ^2U"(p 10 ) = ^3 + 4/pf , a x = + tp^uj 1 / 2 , a\ = 

c< -' 1 ^ 2 yi/2 — ipi/uj 1 / 2 . Note that in the oscillator limit uj\ — > 2, whereas in 
the Wigner limit uj\ approaches the classical result [10]. 

To this order, the spatial wave function is written as e iJ Vo(2/i), where 
■00 ~ ( a i) n |0 >, and 9 is the polar angle associated to the vector p\. Un- 
der a permutation of particles 9 changes by n, thus even \J\ correspond to 
unpolarised spin states, S = 0, and odd \J\ - to polarised states, S = 1. 

The operators hk, k > 3, can be written in general as 

Their matrix elements can be straightforwardly computed. Let us write the 
results for the first nonzero terms of the series (11) in the present case, 

e o = 7P?o-4-> ( 17 ) 
^ P10 

e 2 = u 1 (n + ^), (18) 
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3 1 \ , _ _ 1 



Pio^i \Pio V 4 P 



10 

2 



1 + (30n 2 + 30n + 11), (19) 



Pio^i VPio ' 2 
e 6 = -^^(4/p 4 + l) 4 (2n + l)(47n 2 + 47n + 31) 

oPlO^l 
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, 7 (4/p 4 + l) 2 (6/p 4 + l)(2n + l)(25n 2 + 25n + 19) 



' : (6/p 4 + l) 2 (2n + l)(87n 2 + 87n + 86) 



+ -^L(2n + l)(14n 2 + 14ra + 13) 
Pio^i 

+ ^(2n + l){(10/p 4 + l)(5n 2 + 5n + 9) - 9/p 4 } 
Pio^i 

-— J*— (2n+l), etc. (20) 

n labels the number of excitation quanta. It is a good approximate quan- 
tum number of the problem in the whole interval < (3 < oo. Notice that all 
the 6k with k > 2 go to zero in both the oscillator and Wigner limits. Thus, 
the 1/| J|-expansion provides a nice interpolation scheme. 

In figure la, we compare different approximations to the r.h.s. of (7) 
computed according to the series (11). The excited state with n = 2 and 
| J | = 3 is shown. Convergence is excellent. We can not distinguish among 
the curves containing the corrections 64 and €q (printed as dashed lines) from 
the solid line corresponding to e + 62/] J\- Even at a low value of |J| as 3, 64 
introduces corrections lower than 1%, and eg - lower than 0.3%. The relative 
weight of €q in e = €0 + £2/\J\ + £4/ 'J 2 + e 6/l^| 3 is shown in figure lb. 

In figure 2, the energy \J\e is compared with the estimate obtained in 
[11] from a two-point Pade approximant. We expect this approximant to 
approach the exact energy from above, and to differ from it in no more than 
a few parts in 10 3 [11]. The results are consistent with the expectations. 

Notice that the conditions for the application of perturbation theory are 
fulfilled at any (3 when |J| > 3, i.e. the distances between adjacent levels, 
computed from 60+62/! J|, are much greater than the corrections e±/\ J| 2 + . . .. 
Below, we will see that in the three-body problem the corrections, although 
being small, may cause a partial rearrangement of the spectrum. 
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4 Three electrons 



In the three electron system, there are two distances pi, p 2 , and one angle, 
Q\. The equilibrium configuration is an equilateral triangle (p w = p 2 o, W = 
±7r/2), which side satisfies the equation 

= -l+p? _ 3/ 3 3 pio. (21) 

At its minimum, £/ e //(pio) = 3p 2 /2 — p^~ 2 /2, thus providing a leading 
approximation to the energy 

e o = fp'o " (22) 

Expanding around one of the two equivalent configurations, i.e. fixing 
$io — 7r /2 for example, we obtain a series like (9) in which the first operator 
coefficients are given by 



' d 2 d 2 2 d 2 \ 2i 9 
dy{ dyi p 2 10 dzfj p( dz x 



1/3 
4 Uc 
1 / I 



+ ^ I 1 " — ) ( 5 ^i + 6 ^ 2 + 5 ^ + 3 /4^i ) , (23) 
!6 V Pio/ 



L , 9 a\ 2 . , 9 2 

|si g n(J)(^- y2 ^-i-, 



+ ^-sign(J)(y 2 -y 2 2 )— -— {yl + yl) 

10 

1 A 1 



^J— I 1 r ) ( 19 ?/i + %iV2 + 331/H/2 2 + 9y 3 2 

64pio V Pio/ 

- 9p 2 w y 1 z 2 1 +21p 2 w y 2 z 2 1 ), (24) 
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1/9 9 \ 3 . 2 2 , & 
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^sign(J)( yi 3 - yl)^- + -^(yt + y\) 
Pio 0z i 4 Pio 



1 1 \ /329 4 , 123 9 9 o 9 4 

99 141 41 \ 

-jPwVf z i + 27p? yiy 2 2i 2 + — ploVlzl + -jPw4) , etc ( 25 ) 

The hamiltonian /i 2 , which is to be taken as zeroth order hamiltonian, 



describes a harmonic oscillator of frequency uji = y 3 + 1/ p| in the variable 

2/s = (yi + 1/2) / a/2, plus a combination of harmonic oscillators and a fictitious 
magnetic field in the variables y m = (yi — y 2 )/\/2 and z m = p w Zi/y/2. Note 
that the "magnetic field" comes from the coupling between the angular mo- 
mentum and the variable 9\. We may use a symmetric gauge to describe the 
fictitious field. The hamiltonian is transformed according to h' = e l ^he~ 1 ^ , 
where / = sign(J)y m z m / (2p1 ). h 2 takes the symmetric form 



d 2 1 / 1 \ 2 / 2 10 Id 

K= - 7^ + T 3 + ^ - — + T77T + 




dy 2 4 v pf y ys W £d£ ^ eda 2 

. sign (J) d 



(26) 



where £ = \Jy 2 n + and a = arctan(^ m /?/ m ). The corresponding eigenval- 
ues and eigenfunctions are 



e 2 = uJi(n s + 1/2) + cu 2 (2n + \m\ + 1) + sign( J)mc;3, (27) 
^0 ~ //„ s (vW2y s )e- Wl ^ /4 e |m| e- W252/4 4 m| (^ 2 /2)e ma , (28) 

in which c<j 2 = \Js/2 — pio/2, ^3 = 1/ Pio, an d n an d m are respectively prin- 
cipal and magnetic quantum numbers in the fictitious field. In the oscillator 
limit: uj\ — > 2, uj 2 , ^3 — ► 1, whereas in the Wigner limit: pio — > 00, c<Ji — > \/3, 
cj 2 — > ^3/2, and cj 3 — > 0. At finite /3 there is no degeneracy between the 
states of h 2 , thus we may use nondegenerate perturbation theory to compute 
corrections. 

We shall identify which of the levels (J,n s ,n,m) may describe electron 
states. This set of quantum numbers will serve as approximate quantum 
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numbers for the three-electron problem. At intermediate couplings, a small 
mixing between the states of h' 2 will be induced by the corrections h' 3 and h' 4 . 

Let us recall the symmetry requirements for three-electron wave functions. 
A spatially antisymmetric wave function, = e lJ ~ipo, where S accounts for 
global rotations, corresponds to a spin-polarised state, S = 3/2, whereas 
a mixed-symmetry \[> is related to an unpolarised spin state, S = 1/2. 
Symmetry transformations shall be expressed in terms of the variables we 
are using, i.e. y s , £ and a. Additionally, to identify physical states, \I/o 
shall be compared with oscillator functions at (3 — > 0. A detailed analysis 
is presented in the Appendix of paper [3], to which we refer the reader. 
The main conclusion is that in the (n, m) plane, the spatially symmetric 
and antisymmetric states occupy the lines m = J + 3k , where k is an 
integer. For example, let us consider the first states with J = —3 when 
(3 « 1. The lowest state is (n s , n, m) = (0,0,0). It corresponds to both 
one antisymmetric and one symmetric state. In our scheme, as tunnelling 
effects are not included, they are degenerated in energy. Of course, only 
the antisymmetric state may correspond to a state of three electrons. Very 
near to the lowest state there is the mixed-symmetry doublet, (0,0,1). At 
excitation energies near 2, there are the levels (1,0,0), (1,0,1), (0,0,-1), (0,1,0) 
and (0,1,1). Symmetric and antisymmetric states have m = 0. The next set 
of states is at excitation energies near 4, etc. 

Up to this order, the energy computed from e = e + e 2 /| J\ leads to results 
very similar to that reported in [9]. That is, when f3 ~ 1 the relative error of 
the lowest state with a given \J\ represents only a few percents of the total 
energy at | J| = 3, and decreases considerably as |J| is increased. The reader 
may look at figures 3 and 4 of paper [9] in quality of examples. 

Let us study e when f3 is varied from to oo. As in the two-electron 
problem, the higher corrections, e 4 /\ J\ 2 + . . ., go to zero in both [3 — * and 
(3 — ■> oo limits. However, for certain levels, a rearrangement of the spectrum 
may take place. For example, let us consider spatially antisymmetric states 
with J = —3k, where k is positive and high enough. The lowest state in 
this sector has quantum numbers (0,0,0). We show in figure 3 the excitation 
energies of three states with quantum numbers (0,0,3), (1,0,0) and (0,1,0). 
Their energies are 3(u>2 — 003), 0J1 and 2c^2 respectively. Transitions among 
these states induced by h' 3 , h' 4 , etc. may occur at (3 > 0.6, thus resulting 
in a rearrangement of levels. To compute their energies at the intermediate 
values of f3 where the intersections occur, we shall resolve the degeneracy 
by means of degenerate perturbation theory. On the other hand, there are 
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levels, in particular those with numbers (0,0,0), for which we can apply the 
1/| J|-expansion continuously from /3 — to oo. 

Let us compute the coefficient e 4 for a state with quantum numbers 
(0,0,0). We assume J < and \J\ large enough. These are the relevant 
states at high magnetic fields. We shall compute the matrix elements enter- 
ing the expression 



e 4 = < 0,0,0|/4|0,0,0 > 

_ < 0; 0; 0|m^s; '"- "-S; "-; '"-|"-3l"; » ^ 



< 0, 0, 0j/i^| 


n s , n, m >< n s , n, m\ 


K\ 


0,0,0 > 


n s ui + (2n + 


m\)u)2 — moo?, 



Taking into account the explicit form of h' 3 , i.e. 

d 



where 



ti a = A-jj- + By, + Cyl + D, (30) 



A = (31) 

PlO 



3 / 1 



PlO VOZ m p w OZ m 4p 10 



16 V AJ [vl + zl) ' (32) 



4pio V Pio, 



C = -T-( 1 + ^) ? (33) 



D = -^"f 1 - ^H^-3^4), (3D 
-J^Pio V Pio, 



we obtain that the sum over intermediate states in (29) contains only a few 
terms. The following result for e 4 arises, 



C 2 

e 4 = < o, 0, OKIO, 0, > < 3\y 3 s \0 > 

3ui 



1 , , o/l 1 \ A 2 . d . 
- -<0,3D0,0> 2 + ■ +— <1 — 0> 

3 \oj<2.-uj3 U2 + U3J u)i oy s 
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2C 

- — < l\y s \0 >< l\y 3 s \0 >< 0,015(0,0 > 

C 2 1 

- — < llyflO > 2 < llyJO > 2 < 0,015(0,0 > 2 

„, Ift 2 f < 0,2(5(0,0 > 2 < 0,-2(5(0,0 > 2 1 

- <ly s 0> 2 < ' ' r + J 35 

{ UJi + 2{u 2 - UJ 3 J UJ\ + 2(^2 + W3) J 



where 



< 0,0,0|/i 4 |0,0,0 > 





1 \ 15 














Va;2 










(36) 



< 0,3(5(0,0 > = --^ 1-t , (37) 

< 0,0(5(0,0 > = (38) 



2pio ' 

1 002 1 

Pio 4 Pio 4w 2 Pi : 



<0, ±2(5(0, 0> = ±-3 3-, (39) 



and the matrix elements involving y s are ordinary oscillator matrix elements, 
that is < l\y s \0 >= l/y/uJ[, etc. Note that, as mentioned, e 4 — > in both 
the (3 — > and (3 — > 00 limits. 

In figure 4 we compare the 1/| J (-expansion and the two-point Pade ap- 
proximant P^{f3) of [11] for the lowest antisymmetric state with J = —3. 
As in the two-electron case, the relative error is of the same order of the 
incertitude of the Pade estimate, i.e. a few parts in 10 3 . We show in figure 
5 the relative weight of e 4 in e = eo + ^/\J\ + e 4 /</ 2 - For states with still 
higher values of | J|, the relative magnitude of the contribution of e 4 decreases 
considerably. 



12 



5 Conclusions 



We have computed the energy levels of two and three electrons in an ideal 
parabolic two-dimensional quantum dot in the presence of a magnetic field. 
States with high angular momentum and a low number of excitation (vi- 
brational) quanta were considered. The energy was found as a series in the 
"small" parameter 1/|J|. The series exhibits good convergence properties 
even at | J| = 3. 

At still lower values of the angular momentum the method is not ap- 
plicable. However, we can obtain reliable analytic estimations by means of 
an alternative procedure: the two-point Pade approximants, which are ex- 
tremely simply constructed for the low- lying states [11]. 1/| J|-expansions 
and Pade approximants were shown to give similar results at intermediate 
| J | values, i.e. at \J\ — 3. 

Other problems are currently being treated along these lines. For exam- 
ple, the effects of anyonic statistics on the energy levels of two quasiparticles 
in a quantum dot may be straightforwardly studied by means of the 1/| J|- 
method [12]. 
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Figure Captions 



Figure 1. 

a) The energy of the | J] = 3, n — 2 state of two electrons. 

b) The relative weight of €q in e. 

Figure 2. Relative difference between \J\e and the Pq^ Pade approximant 
for two electrons in a state with | J| = 3, n — 0. 

Figure 3. Excitation energies of three states with J = —3k, k - positive and 
high enough. Quantum numbers are indicated. 

Figure 4. Relative difference between \J\e and the P 5i4 Pade approximant 
for three electrons. The lowest antisymmetric state with J = — 3 is 
studied. 

Figure 5. Relative weight of e 4 in e. The same state as in figure 4 is 
considered. 
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